Astrophysics and Space Science 
DOI 10.1007/sl0509-010-0356-4 



EXOTIME: searching for planets around pulsating 
subdwarf B stars 



o- 

c3 



o\ 



m 

Oh- 
I 

o 



en 

O 

o 



13 



Sonja Schuh 1 ' 2 • Roberto Silvotti 3 ' 4 • Ronny Lutz 2 ' 5 

• Bjorn Loeptien 2 • Elizabeth M. Green 6 • 

Roy H. 0stensen 7 • Silvio Leccia 4 • Seung-Lee Kim 8 

• Gilles Fontaine 9 • Stephane Charpinet 10 • 
Myriam Francoeur 9 • Suzanna Randall 11 • 
Cristina Rodriguez-Lopez 10 ' 12 • 

Valerie van Grootel 10 • Andrew P. Odell 13 • 
Margit Paparo 14 • Zsofia Bognar 14 • 
Peter Papics 7 ' 15 • Thorsten Nagel 1 • 
Benjamin Beeck 2 • Markus Hundertmark 2 • 
Thorsten Stahn 5 • Stefan Dreizler 2 • 
Frederic V. Hessman 2 • Massimo Dall'Ora 4 • 
Dario Mancini 4 • Fausto Cortecchia 4 • 
Serena Benatti 16 • Riccardo Claudi 17 • 
Rimvydas Janulis 18 

• Springer- Verlag 2010 



Sonja Schuh 
Roberto Silvotti 
Ronny Lutz 
Bjorn Loeptien 
Elizabeth M. Green 
Roy H. 0stensen 
Silvio Leccia 
Scung-Lee Kim 
Gilles Fontaine 
Stephane Charpinet 
Myriam Francceur 
Suzanna Randall 
Cristina Rodriguez-Lopez 
Valerie van Grootel 
Andrew P. Odell 
Margit Paparo 
Zsofia Bognar 
Peter Papics 
Thorsten Nagel 
Benjamin Beeck 
Markus Hundertmark 
Thorsten Stahn 
Stefan Dreizler 
Frederic V. Hessman 
Massimo Dall'Ora 
Dario Mancini 
Fausto Cortecchia 
Serena Benatti 
Riccardo Claudi 
Rimvydas Janulis 



x Institut fur Astronomie und Astrophysik, Kepler Center for As- 
tro and Particle Physics, Eberhard-Karls-Universitat, Sand 1, 
72076 Tubingen, Germany 

2 Georg-August-Universitiit Gottingen, Institut fur Astrophysik, 
Friedrich-Hund-Platz 1, 37077 Gottingen, Germany 
3 Istituto Nazionale di Astrofisica (INAF), Osservatorio Astro- 
nomico di Torino, Strada Osservatorio 20, 10025, Pino Torincsc, 
Italy 

4 Istituto Nazionale di Astrofisica (INAF), Osservatorio Astro- 
nomico di Capodimonte, via Moiariello 16, 80131 Napoli, Italy 
5 Max-Planck-Institut fur Sonnensystemforschung, Max-Planck- 
Strae 2, 37191 Katlenburg-Lindau, Germany 

6 Steward Observatory, University of Arizona, 933 North Cherry 
Avenue, Tucson, AZ 85721, USA 

7 Instituut voor Sterrenkunde, K. U. Leuven, Celestijncnlaan 
200D, 3001 Leuven, Belgium 

8 Korca Astronomy and Space Science Institute, Dacjcon 305-348, 
Korea 

9 Departement de Physique, Universite de Montreal, C.P. 6128 
Succ. Centre- Ville, Montreal, H3C 3J7, Canada 
10 Laboratoire d'Astrophysique de Toulouse- Tarbes, Universite 
de Toulouse, CNRS, 14 Avenue Edouard Belin, 31400 Toulouse, 
France 

11 ESO, Karl-Schwarzschild-Str. 2, 85748 Garching bei Miinchen, 
Germany 

12 Departamento de Ffsica Aplicada, Universidade de Vigo, Vigo 
36310, Spain 

"Northern Arizona University, Flagstaff, AZ 86011, USA 
14 Konkoly Observatory of the Hungarian Academy of Sciences, 
P.O. Box 67., H-1525 Budapest XII, Hungary 
15 E6tvos Lorand University Faculty of Science, Pazmany Peter 
setany 1/A, H-1117 Budapest, Hungary 

16 CISAS, Universita degli studi di Padova, via Venezia 15, 35131 
Padova, Italy 

17 Istituto Nazionale di Astrofisica (INAF), Osservatorio Astro- 
nomico di Padova, vicolo Osservatorio 5, 35122 Padova, Italy 
"institute of Theoretical Physics and Astronomy, Vilnius Uni- 
versity, 12 A. Gostauto Street, 01108 Vilnius, Lithuania 



Abstract In 2007, a companion with planetary mass 
was found around the pulsating subdwarf B star 
V391 Pegasi with the timing method, indicating that a 
previously undiscovered population of substellar com- 
panions to apparently single subdwarf B stars might 
exist. Following this serendipitous discovery, the 
EXOTIMBJ monitoring program has been set up to 
follow the pulsations of a number of selected rapidly 
pulsating subdwarf B stars on time-scales of several 
years with two immediate observational goals: 

1) determine P of the pulsational periods P 

2) search for signatures of substellar companions 
in O— C residuals due to periodic light travel time 
variations, which would be tracking the central star's 
companion-induced wobble around the center of mass. 

These sets of data should therefore at the same time: 
on the one hand be useful to provide extra constraints 
for classical asteroseismological exercises from the P 
(comparison with " local" evolutionary models) , and on 
the other hand allow to investigate the preceding evolu- 
tion of a target in terms of possible " binary" evolution 
by extending the otherwise unsuccessful search for com- 
panions to potentially very low masses. 

While timing pulsations may be an observationally 
expensive method to search for companions, it sam- 
ples a different range of orbital parameters, inaccessible 
through orbital photometric effects or the radial veloc- 
ity method: the latter favours massive close-in compan- 
ions, whereas the timing method becomes increasingly 
more sensitive towards wider separations. 

In this paper we report on the status of the on- 
going observations and coherence analysis for two of the 
currently five targets, revealing very well-behaved pul- 
sational characteristics in HS 0444+0458, while show- 
ing HS 0702+6043 to be more complex than previously 
thought. 

Keywords Stars: subdwarfs, oscillations, evolution, 
planetary systems, individual: HS 2201+2610, 

HS 0702+6043, HS 0444+0458. 



Note: This paper contains the content of three contributions 
originally presented as 

• Status report on the EXOTIME program 

• EXOTIME photometric monitoring of HS 0444+0458 

• EXOTIME photometric monitoring of HS 0702+6043. 



1 Introduction 



While P measurements exist for a number of pulsat 



ing white dwarfs ('e.g. iKepler et aZ.Nl99lL l2000l 12 005: 



o P 



Costa. Kepler and WingeJ1999HCosta and Keplerl20 08:| 



Mullallv et al\ [2008J) , similar measurements have only 
been published for one pulsating subdwarf B star so 
far, V391 Pegasi (hereafter HS 2201+2610). 

HS 2201+2610 was first di scovered to be a rapidl y 
pulsating subdwarf B sta r by l0stensen et all ( 2001ah . 
Silvotti et al\ (|2002l l2007h were able to derive P values 
for the two strongest pulsation modes, and found an 
additional pattern in the O— C diagrams that revealed 
the presence of a giant planet in a 3.2 year orbit. Ad- 
ditional slow pulsa tions were subsequently reported by 
Lutz et all l |2009bl) . 

From theoretical considerations, both the rapid p- 
mode as well as the slow g-mode oscillations should be 
very stable in the pulsating subdwarf B stars (stable in 
phase, i.e. coherent over time scales of many years, as 
well as stable in amplitude) due to the identical under- 
lying driving: a k mechanism acting on the Z bump. 
Time dependency of phases and amplitudes can how- 
ever be introduced through nonlinear effects; the result- 
ing mode coupling then starts to invalidate the initial 
simple assumption. It is therefore crucial to identify 
target stars that show a stable behaviour. 

Under the right conditions, this provides the possi- 
bility to obtain meaningful results from long-term pho- 
tometric monitoring of rapidly pulsating sdB stars: It 
allows to search for slow changes in pulsation periods 
due to evolutionary effects, testing the time scales in 
sdB evolutionary models, and may be used as an ad- 
ditional constraint for sdB asteroseismological models. 
It also allows to use the star as a clock that provides a 
regular timing signal, which can be exploited to search 
for small periodic deviations from the mean pulsation 
maxima arrival times. Such periodic deviations can be 
caused by a varying light travel time and may hence 
indicate a wobble of the pulsating subdwarf B star's 
location due to one or more low-mass companions, as 
is the case for HS 2201+2610. This timing method is 
sensitive to planetary masses and large period ranges 
not easily accessible with other methods. 

Given the observational efforts required, the detec- 
tion rate of low-mass companions in wide orbits is 
surprisingly high so far. Besides timing of pulsations 
(|Silvotti et aZ.1 120071 ). timing of the sharp eclipses in 
HW Vi r type systems has been used to detect tertiary 



bodies (Lee et al 



2009l:lQian et aZ.ll2009h . Furthermore, 



http://www.na. astro. it/~silvotti/exotime/ 



Geier et all (J2009J) have even successfully employed the 
radial velocity method to detect a close substellar com- 
panion to HD 149382. 
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Table 1 Overview of EXOTIME targets 



target 



coordinates (equinox 2000.) ttib status 



HS 2201+2610 
HS 0702+6043 
HS 0444+0458 
EC 09582-1137 
PG 1325+101 



22:04:12.0 +26:25:07 

07:07:09.8 +60:38:50 

04:47:18.6 +05:03:35 

10:00:41.8 -11:51:35 

13:27:48.6 +09:54:52 



14.3 collecting data planet candidate published, sini unknown 

14.7 collecting data see 13.11 and 111 
15.2 collecting data see 13.21 

15.0 collecting data 

13.8 collecting data 



As it has been suggested that the influence of sub- 
stellar companions or large planetary companions may 
be decisive for the evolution of subdwarf B stars (see 
Discussion in section [6J , the EXOTIME program has 
been set up to find more such objects. We discuss the 
target selection in section [2] However, as amplitude 
variations h ave repeatedly been reported in sdB pul- 
sators (e.g. iKilkennvl 120101) . it is crucial to carefully 
check the stability of the targets first of all. This con- 
tribution mainly discusses the current status of these 
investigations for two of our objects, HS 0444+0458 and 
HS 0702+6043 (section [3]) . 

HS 0444+0458 was first discovered to pulsate by 
( 2001a). and has been further char- 



0stensen et al. 



acterised by 
were discovered 



Reed et al\ (|2007l ) 



Rapid 

by 



oscillations 



Dreizler et 



HS 0702+6043 
(200 2J), simu l taneo us slow oscillations were reported by 
Schuh et all (|2006l ). The on-going EXOTIME observa- 
tions for H S 0702+604 3 have also prev iously been sum- 
marised bv lLutz et al\ (J2008bl l2009a|) - 



2 Observational Data 

2.1 Target selection 

Similar to the criteria applied when first selecting 
HS 2201+2610 for long-term monitoring, we considered 
two sets of characteristics to choose further targets. 
The first set concerns the pulsations: frequencies, am- 
plitudes and phases are obviously required to be stable 
over long periods of time. Although there are early in- 
dicators if these criteria are not fulfilled, it should be 
clear these properties can only be truly verified, in par- 
ticular for the phases, from long-term monitoring data 
in hindsight (see Section [3]). 



Using target properties partly taken from lReed et al. 



(J20071 ). we selected candidates with few but preferably 
more than one or two frequencies, in order to be able 
to fully resolve the frequency spectrum during short 
individual runs while also retaining the possibility to 
obtain an independent O— C measurement from each 
individual frequency. For the amplitudes we leaned to- 
wards medium values, avoiding, where possible, very 



small amplitudes which would have S/N issues, as well 
as very large ones where one might have to start wor- 
rying about no n-linearity e ffects. 

We also used lReed et alx s listing to consider the am- 
plitude stability criterion give n _ by 



Christensen-Dalsgaard. Kieldsen. and Matteil (J200lh . 
which for oscillations similar to solar-like ones identifies 
the stochastically excited pulsations. While this quan- 
titative distinction for the underlying type of driving 
has, not surprisingly at all, been previously shown in- 
adequate for the opacity- driven subdwarf B pulsators, 
we still found it helpful to restrict our list of objects to 
those with a(A)/(A) < 0.5. We took this as a purely 
observational value that overall should favour the more 
stable pulsators. 

The second set of selection criteria was concerned 
with more practical questions, looking for the brighter 
ones of the available candidates and including preferen- 
tially northern hemisphere targets. Compromising be- 
tween all of these considerations, we selected the objects 
in Tabled] as the initial list of targets to be observed. 

2.2 Observations 

Due to not only the different magnitudes and pulsation 
amplitudes as well as the spatial distribution of the tar- 
gets, but also due to the long-term nature of the project 
and the changing availability of telescopes, the observa- 
tions in the data bases for our target objects come from 
a variety of sites, instruments, and observers. The bulk 
of observations for HS 2201+2610 have been pu b lished 
with extensive observing logs by ISilvotti et all (J2002 , 
2Q07). Further data is still being collected and analysed 
in order to confirm the future evolution in the O— C di- 
agram as predicted from the known orbital motion, as 
well as to search for possible further companions at dif- 
ferent orbital periods. 

The two other targets with the best coverage so far 
(after HS 2201+2610) are HS 0702+6043 and 
HS 0444+0458. For these two targets, the full observing 
logs up to date are provided in Tableg](HS 0444+0458) 
and Table [5] (HS 0702+6043), respectively. The obser- 
vations available for these two targets come from tele- 
scopes ranging from 0.5 m to 3.6 m in aperture, and 
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Fig. 1 TNG data on HS 0444+0458: O-C diagram for 
fl 



Fig. 3 TNG data on HS 0444+0458: amplitude variation 
for fl 
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Fig. 2 TNG data on HS 0444+0458: O-C diagram for 
f2 



Fig. 4 TNG data on HS 0444+0458: amplitude variation 
for f2 



were all obtained with CCD detectors, almost all of 
the time through B filters, and with exposure times 
between 5 s and 40 s. 

We note that the observations for HS 0702+6043 
comprise a very large and very high-quality data set 
obtained by Green et al. in 2007/2008 ("Mt. Bigelow 
data set") primarily for astero seismological pu r poses , 
which is presented in detail bv lFrancceur et all ( 20101) 
elsewhere in this issue. We include and examine it here 
exclusively under the aspect of long-term changes. 



Table 2 The two main frequencies in HS 0444+0458 as 
obtained from the full 2008-2009 data set 



ID Frequency [/iHz] 



Period [s] 



amplitude [mmi] 



fl 7311.7522+0.0008 136.76612+0.00016 7.9+0.3 
f2 5902.527 ±0.003 169.4190 ±0.0008 2.4+0.3 



2.3 Data reduction 

Light curves were extracted from (bias-, dark- and flat- 
field-corrected) CCD frames with aperture photometry 
packages, usually directly by the observers. Relative 
light curves were obtained using a prescribed set of ref- 
erence stars for each target. We then applied an extinc- 
tion correction to the relative light curves using low- 
order polynomials. The timestamps of observations are 
given at mid-exposure, and were barycentrically cor- 
rected as well as corrected for leap seconds. The BJD 
is stripped of the leading 24 and therefore appears as 



Table 3 The two main frequencies in HS 0702+6043 as 
obtained from the 2005 Calar Alto data set 



ID Frequency [juHz] 



Period I 



amplitude [mmi] 



fl 2753.9618+0.0047 
£2 2606.0110 ±0.0241 



363.1132+0.006 
383.7282+0.0035 



27.2+0.2 
5.3±0.2 
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Fig. 5 8 nights of Calar Alto data on HS 0702+6043: O-C Fig. 7 8 nights of Calar Alto data on HS 0702+6043: am- 
diagram for fl plitude variation for fl 
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53388 53390 



Fig. 6 8 nights of Calar Alto data on HS 0702+6043: O-C 
diagram for f2 



Fig. 8 8 nights of Calar Alto data on HS 0702+6043: 
plitude variation for f2 



"Reduced" Julian Date in the figures (BRJD). The nor- 
malised relative intensities were given weights l/a^ cs ac- 
cording to the local residual intensity scatter. The sub- 
sequent analysis o f the weighted light curv es was done 
using Period04 bv lLenz and Bregerl (J2005T ). The errors 
for any frequencies, amplitudes, phases derived for vari- 
ous subsets of t he data were calculated according to t he 
prescription bv Montgomery and O'Donoghud (J1999I ). 



3 Coherence 

3.1 Photometric monitoring of HS0444+0458 

We use the six-night data set from TNG in August 2008, 
which is the longest individual observing run and the 
one with the highest S/N available for HS 0444+0458, 



in order to examine the short-term stability of the two 
main pulsations in this star. 

The frequencies used in this exercise were derived 
from the full data set, the values are as listed in Ta- 
ble [5] Keeping these frequencies fixed, we re-derived 
the best mean amplitude for the overall TNG data set, 
as well as the mean phase. Keeping both the frequen- 
cies and the amplitude fixed, we then turned to the 
individual nights to re-determine the phasings. The dif- 
ferences between the mean phase and the nightly phas- 
ings were converted into time lags, and the errors prop- 
agated from the error of the mean phase and the errors 
of the nightly phasings. The results for the frequencies 
fl and f2 are shown in Fig. [T]and Fig. [5J 

There are no serious discontinuities in the phases of 
the pulsations that would have to be attributed to insta- 
bilities in the pulsational behaviour of HS 0444+0458; 
all deviations are consistent with the uncertainties ex- 
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Fig. 9 60 nights of Mt Bigelow data on HS 0702+6043: 
O — C diagram for fl 



Fig. 11 60 nights of Mt Bigelow data on HS 0702+6043: 
amplitude variation for fl 



-20 t 

54400 




S 0.006 

Q. 

E 



0.003 b 

54400 



54450 54500 

BRJD 



Fig. 10 60 nights of Mt Bigelow data on HS 0702+6043: 
O — C diagram for f2 



Fig. 12 60 nights of Mt Bigelow data on HS 0702+6043: 
amplitude variation for f2 



pected from the measurement uncertainties alone (short 
individual runs equal relatively large uncertainties in 
the phasing compared to phasings derived from an en- 
semble of several nights, as used for the construction of 
long-term O— C diagrams). 

The reasonably flat O— C diagrams therefore show 
that HS 0444+0458's pulsations are coherent on this 
time scale. Extrapolating this result at the given over- 
all frequency precision, the cycle count remains reliable 
during the 2-3 month gaps in between the data sets 
obtained for HS 0444+0458. 

Similarly, Fig. [3] and Fig. @] show the degree of varia- 
tion in the amplitudes for fl and f2 from night to night 
during the TNG run. 



3.2 Photometric monitoring of HS 0702+6043 

The coherence and amplitude stability in HS 0702+6043 
were tested in two independent different data sets. 

The first data set consists of eight nights of data 
obtained at Calar Alto in 2005. The two main rapid 
pulsations we re recovered at the values given in Table [3] 
(compare also lLutz et q/.ll2008a ). They were used in the 
same way as in the previous section to construct O— C 
diagrams (Fig. [5] and Fig. [6]) as well as to investigate 
the amplitude stability (Fig. [7] and Fig. [5} for fl and f2 
in HS 0702+6043. 

Despite the longer data set and higher S/N (due 
to both a brighter target and higher pulsation am- 
plitudes) and hence smaller error bars, the phases in 
HS 0702+6043 appear more scattered than in 
HS 0444+0458, and the amplitude in particular for the 
main frequency fl significantly less stable. 
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The second data set co mprises the core 60 n ights of 



the Mt. Bigelow run (see iFrancceur et aLll2010l for the 



frequency analysis). The most immediate and notable 
result from this data set is that HS 0702+6043 shows 
a lot more frequencies both in the p- as well as in the 
g-mode domain than previously published. 

We treated the data independently as above to ob- 
tain near-continuous O— C diagrams spanning more 
than 100 days for fl and f2 (Fig. [9] and Fig. [TDJl as 
well as diagrams of the evolution of the corresponding 
pulsation amplitudes (Fig. [11] and Fig. [12]) during this 
time. The results for the amplitude var iation confirm 
those obtained by IFrancceur et al\ (|2010l) . 

All of these figures reveal a rather complex behaviour 
that needs to be discussed in detail. The O— C diagram 
for fl shows a significant trend which looks more or less 
linear (which is also seen in the associated amplitude). 
Normally this would be a telltale sign for a mean fre- 
quency chosen at a slightly wrong value, which can be 
corrected to subsequently get rid of the trend. How- 
ever, all attempts to "improve" the frequency have so 
far proven unsuccessful to eliminate this trend (a result 
consistent with the mean solution for the full 2004-2009 
data), so it must represent a higher-order effect. 

The O—C diagram for f2 shows an even more com- 
plex behaviour, with roughly oscillatory shape (with a 
very similar evolution in the amplitude), indicative of 
beating. With less than two full "beating cycles" cov- 
ered even by this very extensive data set, however, the 
hypothetical multiplet structure potentially generating 
this beating still cannot be adequately resolved. 

It is worth noting that the O— C diagrams and am- 
plitude variation results do not change when increas- 
ingly more frequencies are included simultaneously in 
the analysis (up to 18 frequencies tested); this does not 
include, however, hypothetical additional very close fre- 
quencies within the frequency resolution of frequencies 
already considered. 

It would therefore be tempting to claim that this 
data set reveals first indications for what the EXOTIME 
program is, after all, looking for: drifts or even peri- 
odic residuals in the O— C diagrams that can only be 
explained by the influence of substellar companions. 
However, at most one of the trends in the O— C dia- 
grams could be attributed to one or more unseen com- 
panion, as the light travel time would have to change 
in exactly the same way and with exactly the same am- 
plitude for all pulsations, independently of their other 
properties. Such a similarity in the patterns of the O— C 
diagrams for different pulsations has however not been 
uncovered in this analysis so far, and all remaining ir- 
regularities would still have to be explained differently. 
In the meantime, it remains unclear if external effects 



(including unresolved close frequencies) causes these 
features or if the pulsations truly drift very slowly in 
amplitude, phase or period. 

Finally, the fact that the O—C diagrams for the much 
shorter Calar Alto data alone does not look so bad (al- 
though not as good as that for HS 0444+0458) by itself, 
while the O—C diagrams for the much longer continu- 
ous Mt. Bigelow data set proves the star's pulsations to 
be much more complicated, reminds us that all conclu- 
sions derived from shorter runs should be interpreted 
with a fair amount of caution. 

3.3 Follow-up observations of HS 2201+2610 

The O—C diagrams, as well as pulsational amplitudes, 
derived for HS 2201+2610 have been published else- 
where. An extension of the diagrams for fl and f2 of 
this star is work in progress, with significant amounts 
of new data available by now. In addition to the con- 
tinuing long-term photometric monitoring data in one 
filter, multi-colo ur data and t i me-re solved spectra have 
been obtained. ISchuh et all (J2009I ) sketch an idea of 
how to combine such observations to get an estimate 
of the orbital inclination i for the companion, which 
would allow to transform the current rri2 sin i measure- 
ment into a true mass m^. We note that efforts in 
this direction have however only indicated so far that 

a) HS 2201+2610 is probably a slow rotator and that it 

b) has very small pulsational radial velocity amplitudes. 



4 Long-term O — C diagrams 

While the long-term O-C diagrams for HS 2201+2610 
are available elsewhere, their construction for 
HS 0702+6043 is work in progress, and is currently still 
hampered by the difficulties detailed in section [3~2l 

We therefore focus on presenting the preliminary 
O-C diagrams for HS 0444+0458 here. The overall cov- 
erage (a total of 8 month of data) is much shorter here 
in comparison to HS 2201+2610 and HS 0702+6043, so 
that obviously no results for P or even more complex 
residuals can be presented yet. 

The very important preliminary result for this star, 
however, is that the O—C diagrams for both its fre- 
quencies fl and f2 look fairly flat within the error bars 
over the full 8 month period, as can be seen in Fig. [13] 
and Fig. [H This means that HS 0444+0458 is very 
well suited as a EXOTIME target and definitely worth 
further follow-up. 



6 Discussion 
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Fig. 13 All data on HS 0444+0458: O-C diagram for fl 
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Fig. 14 All data on HS 0444+0458: O-C diagram for f2 

5 Summary 

The main result for HS 0444+0458 so far is that the 
O— C diagram for its main pulsation frequency fl is very 
well-behaved. The O— C diagram for f2 is also approxi- 
mately flat, but with much larger errors and hence less 
significance due to the lower S/N. We conclude that the 
degree of coherence of the pulsations in HS 0444+0458 
confirms it as a good target for long-term monitoring. 
The next few years will show if this reveals a measur- 
able P for the star and/or the signature of an orbiting 
(planetary) companion. 

For HS 0702+6043, previously thought to show a 
relatively s imple spectrum in th e p-mode domain, it 



turns out ( Francceur et all 120101 and the analysis in 



this work) that a significant amount of work must first 
be put into developing a very good understanding of 
the pulsational behaviour before these pulsations can 
be used as a tool to derive P or even search for com- 
panions. 



The increasin g number of detections o f low-mass com- 
panions (ISilvotti et a/.ll2007 : Lee et a l. 2009; Qian et al 



20091 iGeier et q/.ll2009t iGeierlbolol) may indicate the 
existence of a previously undiscovered population of 
substellar and planetary bodies around subdwarf B 
stars. This ha s renewed in terest in the question, ini- 
tially raised bv lSokerl (|1998l ). whether planets might be 
relevant for the formation of subdwarf B stars. 

To discuss the possible influence of a substellar com- 
panion on the primary star's evolution on the red gi- 
ant branch, the above systems must be considered in- 



divid ually. In the case of HS 2201+2610 (jSilvotti et al. 



20071) , for instance, the detected planet is too far out to 
have strongly influenced the envelope ejection process. 
On the other hand, when a star forms planets, it often 
forms full planetary systems with several members, so 
that a remaining detectable far-out planet may serve as 
a tracer of previously present closer-in planets, possibly 
destroyed during the RGB phase. In this sense, the de- 
tection of V391 Pegasi b could be interpreted as backing 
up the suggestion that planets can indeed enhance the 
mass loss on the RGB in such a way that th e subsequen t 
formation of a subdwarf B star is favoured ( Sokcr 201 



The d etection of HD 149382 b (Geier et al. 2009; Geie 



20101) may constitute more direct evidence, with the 
actual object having undergone the common envelope 
phase still present today. 

The situation is different for the tertiary bodies 
in the eclipsin g close binary systems HW Vi r and 
HS 0705+6700 (|Lee et al.ll2009HQian et al.ll2009l) . The 
"problem" of "missing" companions that would be 
needed to explain the formation of the primary sub- 
dwarf B star through one of the accepted binary evo- 
lution channels has never existed here, as suitable sec- 
ondary bodies were already known. 

The detections of low-mass objects around appar- 
ently single sdB stars are therefore the more relevant 
ones when it comes to the discussion whether such com- 
panions can be efficient in forming subdwarf B stars. In 
particular, the amount of angular moment transferred 
and hence the modification of the ejection efficiency 
during a common envelope phase are a matter of debate 
among theorists. In alternative scenarios, the existence 
of planets around subdwarf B stars is explained by for- 
mation processes leading to second generation planets. 

The EXOTIME program has the potential to add 
to the empirical data available on which such a debate 
must ultimately be based. 

Acknowledgements The authors thankfully ac- 
knowledge the contributions by J. Dittmann, A.G. Fay, 
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Table 4 


Photometric data archive for HS0444+0458 










Year 


Month 


Day 


L[h] 


N 


f^rcs 


Exp [s] 


Site a 


Inst. 


Filter 


Observer 15 


2008 


Aug 


19 


0.5 


107 


0.0069 


5 


TNG 3.6m 


DOLORES 


B 


RS, SL 






21 


1.0 


237 


0.0076 


5 


TNG 3.6m 


DOLORES 


B 


RS, SL 






24 


0.7 


173 


0.0060 


5 


TNG 3.6m 


DOLORES 


B 


RS, SL 






26 


1.0 


224 


0.0051 


5 


TNG 3.6m 


DOLORES 


B 


RS, SL 






27 


2.1 


357 


0.0040 


10 


TNG 3.6m 


DOLORES 


B 


RS, SL 






28 


1.7 


230 


0.0042 


10 


TNG 3.6m 


DOLORES 


B 


RS, SL 




Sep 


29 


2.2 


171 


0.0089 


10 


CA 2.2m 


CAFOS 


B 


service (FH) 




Oct 


24 


2.8 


300 


0.0191 


12 


LOAO 1.0m 


2k CCD 


B 


SLK 






26 


4.9 


588 


0.0232 


12 


LOAO 1.0m 


2k CCD 


B 


SLK 






29 


1.2 


163 


0.0175 


12 


LOAO 1.0m 


2k CCD 


B 


SLK 




Nov 


21 


3.2 


263 


0.0092 


10 


CA 2.2m 


CAFOS 


B 


service (MA) 


2009 


Jan 


29 


2.6 


275 


0.0239 


20 


LOAO 1.0m 


2k CCD 


B 


SLK 






31 


4.7 


497 


0.0275 


20 


LOAO 1.0m 


2k CCD 


B 


SLK 




Mar 


20 


1.0 


94 


0.0090 


10 


CA 2.2m 


CAFOS 


B 


service (MA) 






21 


1.3 


135 


0.0121 


10 


CA 2.2m 


CAFOS 


B 


service (MA) 






22 


1.2 


124 


0.0099 


10 


CA 2.2m 


CAFOS 


B 


service (MA) 




Oct 


16 


1.3 


195 


0.0138 


20 


M/N 1.2m 


Ik CCD 


B 


RL 






18 


1.6 


242 


0.0148 


20 


M/N 1.2m 


Ik CCD 


B 


RL 






22 


2.1 


235 


0.0102 


10 


CA 2.2m 


CAFOS 


B 


service (AA) 






23 


2.9 


330 


0.0082 


10 


CA 2.2m 


CAFOS 


B 


service (MA) 






24 


2.2 


244 


0.0075 


10 


CA 2.2m 


CAFOS 


B 


service (MA) 






25 


1.7 


104 


0.0061 


10 


CA 2.2m 


CAFOS 


B 


service (MA) 






26 


3.4 


387 


0.0086 


10 


CA 2.2m 


CAFOS 


B 


service (MA) 



a CA: Calar Alto Observatory, LOAO: Mt. Lemmon Optical Astronomy Observatory, M/N: MONET/North Telescope, TNG: Tele- 

scopio Nazionale Galileo 

6 FH: Felipe Hoyo, MA: Manuel Alises, RL: Ronny Lutz, RS: Roberto Silvotti, SLK: Seung-Lee Kim, SL: Silvio Lcccia 



J. Laird, CM. Limbach, J. Portouw, M. Revelle, J.M. 
Sierchio, S.M. Story, C. Stratton, A. Strom, P. Wrob- 
lewski (all Steward Observatory), who helped obtain 
the Mt. Bigclow data set. S.S. acknowledges fund- 
ing through the Eberhard-Karls-Universitat Tubingen's 
Teaching Equality program, and also thanks the DAAD 
for allocating a travel grant to attend the Forth Meet- 
ing on Hot Subdwarfs and Related Objects where this 
paper was presented. The work of R.L. on this project 
is funded by a stipend from the IMPRS PhD program 
at MPS Katlenburg-Lindau. C.R.L. acknowledges an 
Angeles Alvarino contract of the regional government 
Xunta de Galicia. Amongst others, the work in this 
paper is 

- based on observations collected at the Ccntro As- 
tronomico Hispano Aleman (CAHA) at Calar Alto, 
operated jointly by the Max-Planck Institut fur As- 
tronomie and the Instituto de Astrofisica de An- 
dalucia (CSIC). 

- based on observations made with the Italian Telesco- 
pio Nazionale Galileo (TNG) operated on the island 
of La Palma by the Fundacin Galileo Galilei of the 
INAF (Istituto Nazionale di Astrofisica) at the Span- 
ish Observatorio del Roque de los Muchachos of the 
Instituto de Astrofisica de Canarias. 



based on data obtained with the MOnitoring NEt- 
work of Telescopes (MONET), funded by the 'As- 
tronomie & Internet' program of the Alfried Krupp 
von Bohlen und Halbach Foundation, Essen, and op- 
erated by the Georg-August-Universitat Gottingen, 
the McDonald Observatory of the University of Texas 
at Austin, and the South African Astronomical Ob- 
servatory. 
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Table 5 


Photometric data archive for HS 0702+6043, part 


a 








Year 


Month 


Day 


L[h] 


N 


•^"rcs 


Exp [s] 


Site a 


Inst. 


Filter 


Observer 5 


1999 


Dec 


06 


5.0 


943 


0.0103 


10 


CA 1.2m 


CCD 


- 


SD, SS 






07 


1.9 


414 


0.0141 


8 


CA 1.2m 


CCD 


- 


SD, SS 






08 


1.4 


241 


0.0085 


10 


CA 1.2m 


CCD 


- 


SD, SS 


2000 


Oct 


08 


0.7 


248 


0.0027 


10 


NOT 2.56m 


ALFOSC 


3W 


R0 


2004 


Feb 


04 


7.3 


685 


0.0233 


30 


Tue 0.8m 


ST7 


G 


TN 






09 


5.8 


244 


0.0016 


60 


StB 2.2m 


2k CCD 


F555W 


EMG 






10 


6.2 


264 


0.0014 


60 


StB 2.2m 


2k CCD 


F555W 


EMG 


2005 


Jan 


17-27 


56.0 


5981 


0.0082 


20, 30 


CA 2.2m 


CAFOS 


B 


SS, TS 


2007 


Dec 


18 


5.9 


498 


0.0442 


40 


Tue 0.8m 


SBIG ST-L 


B 


TN 






19 


6.7 


557 


0.0589 


40 


Tue 0.8m 


SBIG ST-L 


B 


TN 






20 


10.0 


842 


0.0994 


40 


Tue 0.8m 


SBIG ST-L 


B 


TN 






21 


9.3 


760 


0.1350 


40 


Tue 0.8m 


SBIG ST-L 


B 


TN 








12.1 


621 


0.0200 




Goe 0.5m 


SBIG ST-L 


B 


SS, RL, RK 


2007 


Nov 


























424.0 


31638 


0.0027 


50 


MtB 1.55m 


MONT 4k 


S 8612 


EMG & coll. 


2008 


Mar 




















2008 


Feb 


08 


7.2 


601 


0.0124 


40 


Tue 0.8m 


SBIG ST-L 


B 


TN 






09 


7.1 


590 


0.0121 


40 


Tue 0.8m 


SBIG ST-L 


B 


TN 








8.9 


642 


0.0381 


35 


Goe 0.5m 


SBIG ST-L 


B 


RL, BB 






10 


5.8 


488 


0.0134 


40 


Tue 0.8m 


SBIG ST-L 


B 


TN 








12.2 


910 


0.0324 


35 


Goe 0.5m 


SBIG ST-L 


B 


RL, MH 






11 


4.5 


338 


0.0439 


35 


Goe 0.5m 


SBIG ST-L 


B 


RL, MH 






13 


1.1 


86 


0.0749 


35 


Goe 0.5m 


SBIG ST-L 


B 


RL 






15 


5.4 


379 


0.0552 


35 


Goe 0.5m 


SBIG ST-L 


B 


RL, SW 






29 


3.7 


360 




20 


Loi 1.5m 




B 


service (ADB) 




Mar 


01 


5.7 


451 




25 


Loi 1.5m 




B 


service (ADB) 






03 


5.1 


429 




25 


Loi 1.5m 




B 


service (ADB) 






05 


0.6 


47 


0.0238 


35 


Goe 0.5m 


SBIG ST-L 


B 


RL, SS 






13 


3.7 


126 


0.0069 


30 


Asi 1.8 


Afosc 


R 


SB, RC 








2.3 


248 


0.0105 


30 


Kon 1.0m 


PI VA 1300B 


B 


MP, PP 






15 


4.6 


597 


0.0095 


25 


Kon 1.0m 


PI VA 1300B 


B 


MP, PP 






18 


3.7 


348 


0.0111 


35 


Kon 1.0m 


PI VA 1300B 


B 


MP, PP 






19 


1.1 


137 


0.0088 


25 


Kon 1.0m 


PI VA 1300B 


B 


MP, PP 




May 


11 


3.0 


204 


0.0315 


40 


Goe 0.5m 


SBIG ST-L 


B 


SW 






12 


1.2 


160 


0.0380 


40 


Goe 0.5m 


SBIG ST-L 


B 


SW 






13 


3.2 


266 


0.0254 


40 


Tue 0.8m 


SBIG ST-L 


B 


AH 




Oct 


15 


1.3 


113 


0.0085 


10 


CA 2.2m 


CAFOS 


B 


service (LM) 






16 


0.8 


65 


0.0101 


10 


CA 2.2m 


CAFOS 


B 


service (MA) 






17 


0.5 


48 


0.0152 


10 


CA 2.2m 


CAFOS 


B 


service (MA) 






18 


1.6 


143 


0.0068 


10 


CA 2.2m 


CAFOS 


B 


service (MA) 






25 


2.6 


285 


0.0089 


30 


Kon 1.0m 


PI VA 1300B 


B 


ZB 








4.4 


518 


0.0107 


15 


LOAO 1.0m 


2k CCD 


B 


SLK 






28 


2.9 


326 


0.0104 


15 


LOAO 1.0m 


2k CCD 


B 


SLK 






30 


3.9 


455 


0.0104 


15 


LOAO 1.0m 


2k CCD 


B 


SLK 




Nov 


21 


3.0 


254 


0.0048 


10 


CA 2.2m 


CAFOS 


B 


service (MA) 






23 


5.7 


496 


0.0064 


10 


CA 2.2m 


CAFOS 


B 


service (MA) 




Dec 


16 


3.0 


243 


0.0068 


25 


M/N 1.2m 


Ik CCD 


B 


RL, BL 






20 


2.7 


202 


0.0116 


25 


M/N 1.2m 


Ik CCD 


B 


RL 


2009 


Jan 


06 


1.6 


143 


0.0364 


30 


Goe 0.5m 


SBIG ST-L 


B 


RL, US 






22 


3.7 


334 


0.0080 


25 


M/N 1.2m 


Ik CCD 


B 


RL 






28 


5.4 


473 


0.0147 


20 


LOAO 1.0m 


2k CCD 


B 


SLK 






30 


6.6 


662 


0.0138 


20 


LOAO 1.0m 


2k CCD 


B 


SLK 



"see footnote in Table [6] for the site abbreviations 

6 see footnote in Table [6] for the observer abbreviations 



EXOTIME: searching for planets around pulsating subdwarf B stars 



11 



Table 6 Photometric data archive for HS 0702+6043, part b 



Filter Observer 15 



Year Month Day L [h] 



N 



Exp [s] Site a 



Inst. 



2009 



Feb 



Mar 



Sep 



Oct 



18 


2.4 


202 


0.0256 


30 


23 


3.6 


304 


0.0126 


25 


25 


6.4 


564 


0.0076 


25 


26 


6.5 


536 


0.0111 


25 


27 


6.6 


556 


0.0083 


25 




5.7 


532 


0.0067 


15 


01 


6.2 


497 


0.0090 


25 


18 


3.3 


279 


0.0060 


25 


19 


3.6 


291 


0.0098 


25 


20 


2.4 


243 


0.0058 


10 


21 


2.4 


232 


0.0068 


10 




6.0 


351 


0.0150 


30 


22 


2.4 


245 


0.0049 


10 




4.9 


279 


0.0140 


30 


24 


3.4 


243 


0.0092 


25 


26 


2.8 


251 


0.0081 


25 




3.9 


396 


0.0118 


20 


27 


1.3 


160 


0.0185 


20 


28 


3.9 


404 


0.0172 


20 


29 


2.9 


309 


0.0112 


20 


30 


3.8 


386 


0.0144 


20 


31 


4.0 


399 


0.0154 


20 


27 


3.9 


406 


0.0132 


25 


28 


4.0 


420 


0.0117 


25 


16 


4.0 


412 


0.0099 


25 


18 


0.9 


100 


0.0109 


25 


19 


6.1 


645 


0.0119 


25 


22 


3.9 


359 


0.0076 


10 


23 


3.6 


334 


0.0067 


10 


24 


3.1 


289 


0.0061 


10 


25 


2.8 


297 


0.0082 


25 




3.4 


315 


0.0078 


10 


26 


3.0 


279 


0.0064 


10 



Goe 0.5m 


SBIG ST-L 


B 


RL, TOH 


M/N 1.2m 


Ik CCD 


B 


RL, US, MH 


M/N 1.2m 


Ik CCD 


B 


RL 


M/N 1.2m 


Ik CCD 


B 


RL, LN 


M/N 1.2m 


Ik CCD 


B 


RL 


Asi 1.8m 


Afosc 


no 


SB, RC 


M/N 1.2m 


Ik CCD 


B 


RL 


M/N 1.2m 


Ik CCD 


B 


RL 


M/N 1.2m 


Ik CCD 


B 


RL 


CA 2.2m 


CAFOS 


B 


service (MA) 


CA 2.2m 


CAFOS 


B 


service (MA) 


Asi 1.8m 


Afosc 


B 


SB, RC 


CA 2.2m 


CAFOS 


B 


service (MA) 


Asi 1.8m 


Afosc 


B 


SB, RC 


M/N 1.2m 


Ik CCD 


B 


RL 


M/N 1.2m 


Ik CCD 


B 


RL, US 


LOAO 1.0m 


2k CCD 


B 


SLK 


LOAO 1.0m 


2k CCD 


B 


SLK 


LOAO 1.0m 


2k CCD 


B 


SLK 


LOAO 1.0m 


2k CCD 


B 


SLK 


LOAO 1.0m 


2k CCD 


B 


SLK 


LOAO 1.0m 


2k CCD 


B 


SLK 


M/N 1.2m 


Ik CCD 


B 


RL 


M/N 1.2m 


Ik CCD 


B 


RL 


M/N 1.2m 


Ik CCD 


B 


RL 


M/N 1.2m 


Ik CCD 


B 


RL 


M/N 1.2m 


Ik CCD 


B 


RL 


CA 2.2m 


CAFOS 


B 


service (AA) 


CA 2.2m 


CAFOS 


B 


service (MA) 


CA 2.2m 


CAFOS 


B 


service (MA) 


M/N 1.2m 


Ik CCD 


B 


RL 


CA 2.2m 


CAFOS 


B 


service (MA) 


CA 2.2m 


CAFOS 


B 


service (MA) 



"Asi: Asiago 182cm Copcrnico Telescope, CA: Calar Alto Observatory, Goe: Gottingcn IAG 50cm Telescope, Kon: Konkoly RCC 
Telescope, LOAO: Mt. Lcmmon Optical Astronomy Observatory, Loi: Loiano 152cm Telescope, M/N: MONET/North Telescope, 
MtB: Mt. Bigclow Kuiper Telescope, NOT: Nordic Optical Telescope, StB: Steward Observatory Bok Telescope, Tue: Tubingen 80cm 
Telescope 
AA: Alberto Aguirre, ADB: Antonio Dc Blasi, AH: Agnes Hoffmann, BB: Benjamin Beeck, BL: Bjorn Loeptien, EMG: Elizabeth 
M. Green, LM: Luzma Montoya, LN: Lisa Nortmann, MA: Manuel Alises, MH: Markus Hundcrtmark, MP: Margit Paparo, PP: Peter 
Papics, RC: Riccardo Claudi, RK: Rcnatc Kruspc, RL: Ronny Lutz, R0: Roy Ostenscn, SB: Serena Bcnatti, SD: Stefan Drcizlcr, 
SLK: Seung-Lee Kim, SS: Sonja Schuh, SW: Sascha Werhahn, TN: Thorsten Nagel, TOH: Tim-Oliver Husser, TS: Thorsten Stahn, 
US: Ulf Seemann, ZB: Zsofia Bognar 
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